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Abstract 



It is shown that large perturbative corrections found previously for semileptonic 
beauty and charm decays are associated with using inappropriate pole masses. The 
latter, in the perturbative expansion, suffer from the I/tuq infrared renormalon 
which is absent in the widths, which leads to similar large corrections in mg. Pole 
masses are neither measured directly in experiment. If the widths are related to 
parameters determined in experiment, the overall impact of the calculated second 
order corrections becomes strongly suppressed and leads to less than 1% change in 
\Vcb\ and \Vub\- Even in charm decays the perturbative corrections appear to be 
very moderate in the consistent OPE-compliant treatment. The updated estimate 
of \Vcb\ is given, based on recent accurate determination of rub and ^^(l GeV). The 
theoretical accuracy of determination of |\4b| from rsi(6 u) appears to be good 
as well. 



* Permanent address 



Inclusive semileptonic decays of B mesons provide at present the most accurate 
determination of the CKM mixing element \Vcb\ and, if one could measure the inclu- 
sive width r(i? — s> X„£z/), the value of \Vub\ as well. In order to reach the ultimate 
possible theoretical precision it is practically necessary to account for purely per- 
turbative corrections as accurately as possible. The first order corrections are taken 
directly from the QED calculations for fi decay known since 50s [|l|. Recently the 
part of the second order corrections has been computed, which is associated with 
the running of the strong coupling (we will refer to such approximation below as to 
the BLM approach). Let us introduce the following general notations: 



r°i(mg/m6) \ vr V vr 



Then the result for rsi(6 qiu) reads (for as{mf,) in the V scheme and 
Uf = 3) 

ai = —2.41 , 02 = —19.7 at mq/nib = 

Oi = —1.67 , 02 = —8.8 at rrig/mh = 0.3 (2) 

where the first line corresponds to the h ^ u decays and the second one is relevant 
for the b ^ c transitions. In Eq. (^ only the BLM part of 02 is shown. In the 
absence of the complete calculations we will discuss only this part to the rest of this 
paper without explicit reminding this reservation. It is worth noting that the values 
of a2 quoted in Eq. @ correspond to the V scheme which has the direct physical 
meaning in the BLM calculations For this reason we will consistently use below 
the V scheme. 

The above values of 02 are surprisingly large at first sight and could be thought to 
change significantly the theoretical estimates for the semileptonic widths, in partic- 
ular in the 6 — > u decay, once (a^/Tr)^ terms are accounted for. It appears, however, 
that the actual impact of these terms is essentially smaller, because similar large 
(a^/vr)^ corrections affect the values of the heavy quark masses. This "conspiracy" 
is the reflection of the dominance of the leading l/mg infrared renormalon which, 
in reality, is absent in the widths once they are expressed in terms of observable 
quantities [Q. The numerical situation with perturbative corrections in semilep- 
tonic beauty decays is, therefore, rather similar to the case of t —>■ b + W width 



1 Heavy Quark Masses and Perturbative Correc- 
tions to Widths 

Heavy quark masses are not directly observed and must be determined indirectly 
from independent measurements. At present the most accurate determination of 

^MS scheme is unnatural here; in this scheme one would have aj^^"^^ = a^2^^ + 56/12 ai where 
6=11 — 2/3n/~9is the first coefficient in the strong coupling /3-function. 
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rrih follows from the sum rules for bb threshold production in e+e^ annihilation [Q; 
the value of nih there undergoes perturbative corrections as well and this must be 
accounted for properly. 

The typical momentum scale in the analysis of moments of the spectral density 
in Ref. P is about 2 GeV, and, therefore, these sum rules are sensitive to the mass 
normalized at the scale /x > 1 GeV. It does not prevent expressing the final result 
in terms of the "one loop pole mass" 



(1) / N —dmiu) , , OLsilA /o\ 
"^poie - ^KN 3 ^^ - + Cm A* (3) 

^ dfl TT 



as long as the one loop corrections are concerned. However, accounting for the 
terms shifts the corresponding value of the pole mass to the value of the "two loop 
pole mass" which is given in the BLM approximation by 



TT 



Numerically this increase is as large as 150 MeV and this strongly affects the width 
as well. (Eqs. (H) and (§) assume that /i <C mq.) 

It is more practical, both theoretically and phenomenologically, to reexpress 
even the perturbative relation for the width in terms of the running masses normal- 
ized at the scale /i. To the second order in it would require the explicit knowledge 
of the complete two loop scale dependence of the heavy quark mass, which is not 
known yet 0. Being interested only in the BLM-type contribution to a2, on the other 
hand, we need to account only for the terms ~ b( ""^^^ y which are easily obtained 
from the first order calculation. This relation becomes especially simple in the heavy 
quark expansion which holds when /i -C mq : 

mpoic ^ rn{fi) + Crn'^^^^ 1^ ^log ^ + 1^ yU + ... (5) 

where A is an arbitrary renormalization scale for as (in the V scheme). Then one 
finds for the b ^ u case 



ai(yu) = -- 71" T I + 5 — c 



3 V 4 / rrib 



m 



a2(/i) = aa + 5— c^^ (log- + 1 ) ; (6) 



^The two loop running of mg is known in the dimensional rcgularization, which, however, is 
irrelevant for /i <C rng. mQ{fi) in naive dimensional rcgularization has nothing common to the 
properly infrared defined mass needed for performing OPE B. 
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a(/i) refer to the perturbative coefficients utilizing running masses (for simplicity we 
do not change the argument of as here). Adopting for the purpose of illustration 
/i = 1 GeV and Cm = ^ (see [Q), we get for \ = nit, 

aiin) ^ -1.03 , a2{fi) ~ -3.7 . (7) 

The perturbative corrections appear to be of "normal" magnitude. 

A similar consideration can be applied to the b c decays, with only minor 
technical complications. Here one is relatively close to the SV limit where the 
total width would depend only on mf, — rUc rather than on the absolute values of 
masses. Then the impact of using the pole masses instead of the running ones with 
<ti rric , rrib is suppressed, and one gets essentially smaller values of 02 from the 
very beginning. In particular, if one uses the strong coupling normalized at the scale 
y/mctTif, which is appropriate for the SV kinematics ^ ^ then the value of 02 in the 
V scheme is 

02 ~ -4.3 . 

This coefficient, still somewhat enhanced, per se would lead to the increase in the 
value of \Vcb\ by 2 percentage points. On the other hand, according to Ref. 0, just 
the increase in the input value of mi, by 150 MeV leads to the decrease of \Vcb\ by 
2.5 percentage points. Therefore when one consistently proceeds from the ffist order 
approximation to the second one, the two effects tend to strongly offset each other. 
We will show below that the quoted value of 02 is indeed largely due to the residual 
dependence of the tree level width on the absolute values of masses for actual nic 
and rrih and, therefore, in reality the impact of the second order corrections is further 
suppressed. 

To proceed from the pole masses to mc(/i) and mi,{ii) one can follow the same 
way as for b ^ u, however in this case the leading in ix/rric approximation used 
in Eqs. ®-(0) a priori could be not well justified. The corrections can be readily 
taken into account and, as a matter of fact, appear to be small: the deviation from 
the linear dependence is less than 7% even for n/rric = 1. Let us briefly discuss 
the general case of arbitrary fi/mg. 

In the framework of the BLM approach the exact /i dependence is given by the 
value of the one loop diagram evaluated with the running aglk"^): 

mQ(/i) - mQ(/i') = / FUk') dk' (8) 

where Fm{k'^) is obtained by integrating the expression for the one loop Feynman 
graph over the directions of the gluon momentum in the 4 dimensional Euclidean 
space. In fact, this quantity is conveniently expressed in terms of the one loop cor- 
rection evaluated with the massive gluon propagator: assigning the gluon mass v one 
calculates, say, 5'w!q{v^) substituting 1/k'^ l/(/c^ + i/'^) in the gluon propagator. 
Then k'^F{k'^) for any observable is equal to the discontinuity of the corresponding 
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z/^-dependent quantity with respect to at the negative value of the gluon mass 
squared. For example, for the heavy quark mass one has 



71 



27ii 



-k^ + le) - m^nHi^^ 



-k' - le 



(9) 



Eq. thus expresses the generic function w considered recently in Ref. 110| in terms 
of the gluon mass dependent observable which had been introduced in Refs. 0, [l^, 

Expanding Eq. (||) in as(A) one gets 



"^pole = mQ{^) + 



71 



(10) 

(the pole mass to any finite order is defined as lim^^o mq^fi) ). 

At < rriQ one has |T|, g m^Q\u^) ~ m''Q\o) ~ -fa^ Vu^ and, thus, Fmik"^) ~ 
2/(3v^) which yields Eqs. (^, (^) with Cm = 4/3. The correction factors to these 
approximate relations are easily calculated using the exact one loop expression for 



(1)/ 2\ /n\ 2 «s 

m'Q'{u ) = mQ{0) - 3— "^Q 



^ cia(l + a) log fl + 



1 — a 



The integral is computed straightforwardly (see, e.g. [jl3l) but is inessential for us. 
The discontinuity over u"^ is obtained directly from Eq. (pA]) and takes the form 



1 ^™(t) 1 



2mQ 3m 



^ , , 4 t t3/2 
JF^(t) ~ + 
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(12) 



m: 



Q 



This coincides with the expressions quoted in Refs. |T^, |I0|, |I5[. The correction 
factors corresponding to integration of the exact function J-'^ instead of its approx- 
imate value 4/3 for small fc^, differ from unity by only 0.075 and 0.05 for the first 
and second orders in a^, respectively, even at /.i = rriQ. 
Taking all these corrections into account and using 



ml\Vgb\^ frif 



a 



^0 —1 1 + oi h 02 — + ••• 



TT 



a. 



TT 



(13) 



with 
one gets 
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zo{x) = 1 — 8x — 12x logx + 8x' 
ai(/i) = —1.03 , a2(/i) = —3.7 for b 



u 
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a^(^) = -1.07 , a2(/i) = 0.9 for b^c (14) 

for /X = 1 GeV, rrif, = 4.8 GeV and rric/mb = 0.3 (as previously, I use the V scheme 
as normahzed at rub for b —>■ u and at y/mjrib for b ^ c). The values of the second 
order coefficients are moderate and di{fi) are suppressed as compared to the case of 
the pole masses. 

Let us emphasize that the estimates above were for illustration purposes only; 
to show the dominant source of the perturbative coefficients I deliberately limited 
myself only to the redefinition of masses while not changing the scale of or the 
numerical values of representative quark masses. These additional modifications are 
accounted for properly in the final numerical evaluation; not surprisingly, they do 
not change the width by any noticeable amount. 

We now dwell on the definition of the renormalization point fi for masses used 
above. The separation of low and high momenta can be accomplished technically 
in different ways leading, e.g., to somewhat different values of the coefficient Cm- 
The most natural choice for nonrelativistic expansion is the cutoff over the spacelike 



gluon momenta, \k\ > /x, used in Ref. ||T6[ for heavy flavor transitions in the SV 



kinematics. It is also most natural for the nonrelativistic bb in the T system. This 
method applied to calculation of the heavy quark mass yields the value Cm = 4/3 
used above. In general, in the BLM-type calculations the step-like cutoff over Eu- 
clidean four momentum 6{k'^ ~/^^) acts differently; however, in the case of the heavy 
quark mass /i-dependence the two approaches are equivalent in the leading in /i/mg 
approximation 

It was found in Ref. P|, in accord with the general arguments above, that the one 
loop value of the b quark pole mass mp^jg was strongly correlated with the obtained 
value of as- On the other hand, the one loop pole mass normalized at 1.3 GeV, 

was essentially uncorrelated with as and literally had significantly smaller range of 
variation ||^. Therefore, such a normalization scale can be considered as physically 
appropriate one; at least, it eliminates the strong correlation to the first loop ap- 
proximation. The possible remaining (not enhanced) second order corrections to 
the mass normalized at this scale are then the part of the theoretical uncertainty of 
determination of the heavy quark mass in the approach of Ref. [0. Clearly, one can, 
in principle, use arbitrary normalization point n of approximately this scale. The 
dependence of ai and d2 on /i at = 4/3 is illustrated in Table 1. Perturbative 
corrections appear to be suppressed for any reasonable choice of fi. 

word of caution is worthwhile. If one uses not a Lorcntz invariant cutoff then the relativistic 
dispersion law can be modified. For example, the effective heavy quark mass which determines the 
rest frame energy may start to differ from the one in the dispersion law E(q) — E(0) — /2mQ + ... . 
It is for this reason the slightly different value of Cm = 16/9 was obtained in Refs. |17| T6| considering 
the recoil effects in the second sum rules in the SV limit. I am grateful to M. Voloshin for the 
discussion of this point. In the concrete version of the BLM approach considered in the present 
paper this subtlety is absent. 
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The following remark is in order here. The value of the first order coefficient is 
also noticeably suppressed if one uses the running masses mqi^fj,) above instead of the 
pole masses. Accordingly, the value of 0,2 is less sensitive to the particular scheme 
for as- For example, even in the MS scheme the values of 0,2 for /i = 1.3 GeV become 
—3.4 and —1.5 for b —>■ u and b —>■ c, respectively, i.e. small as well. Moreover, it 
is clear that the value of the whole perturbative factor through the second order, 
1 + diag/n + 02(0^/71)^ does not exhibit any noticeable dependence on the scheme 
chosen for (which is always present at some level due to the truncation of the 
series) . 

The above analysis for the 6 — >■ m case shows that practically the whole large 
second order contribution obtained in Ref. numerically comes from the domain 
of gluon momenta ~ (0.2 0.25)mfe and, therefore, is mainly associated with the 
leading l/mg renormalon in the pole mass which, however, is absent in the width in 
the proper treatment |T^. Applying this observation to the semileptonic charm 
decays we identify that the large "nonconvergent" second order correction pointed 
out in § came, in fact, from the gluon momenta about and below the infrared pole 
in the gluon propagator; this clearly explains the observed numerical behavior. This 
contribution is, therefore, phenomenologically irrelevant and is explicitly excluded 
even from perturbative corrections in the consistent treatment. In particular, for 
~ 400 MeV , which seems to be the minimal reasonable cutoff, the perturbative 
corrections for charm literally become only 



1 - 0.5— - 0.4 ( — ) ~ 1 - 0.07 - 0.008 ~ 0.92 ; 

the second order yields a tiny contribution only. Altogether one has very moder- 
ate effect ip. The non-BLM higher order terms are potentially more important, 
but most likely are by far dominated by calculable power nonperturbative [13] and 



2 



unknown "exponential" terms |I9 



We see that in the conventional approach, when all intermediate results are 
expressed in terms of the pole mass, the final result for the width is rather sensitive 
to the treatment of the infrared region in the Feynman integrals; already at an 
accuracy level of 5 10% in beauty decays care must be taken to do it in a self- 
consistent way. On the other hand, just the physics at this scale affects the widths 



very little [|T8], ^ (as long as one keeps the high scale heavy quark masses fixed 
rather than the masses of observable hadrons). Therefore it is advantageous to get 
rid of this low energy region and, therefore, of the pole mass in the extraction of 
weak mixing parameters from start. It is very important from phenomenological 
point of view that the sum rule analysis of the bb production allows the direct 
determination of the high scale mb{fi) with in the interval ~ (0.75 -t- 1.5) GeV. 
The theoretical accuracy there is optimal for /i ^ 1 GeV, in accord with the general 
OPE picture, where the impact of low momentum region is described by the vacuum 
expectation value of ^^^(/i) and was estimated to be negligible for any reasonable 
normalization scale. 
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This fact is crucial for rehabihty of theoretical calculations, because the per- 
turbative treatment in the region where is not small is rather involved and, as 
a matter of fact, is sensitive to the various approximations. For example, going 
beyond the BLM computations can change drastically all corrections coming from 
this domain, both for the pole mass and for the width expressed through it. The 
situation is quite different if the low momentum region is explicitly excluded. In 
this case the higher order terms produce small impact. This will be numerically 
illustrated in Sect. 2 when the result for \Vqb\ changes only by a minute amount due 
to the second order corrections, when expressed in terms of ^q(1 GeV). Only then 
one can count on smallness of higher order corrections, and on the moderate size of 
corrections due to the residual terms not captured in the BLM approximation. 

We need to emphasize here two points related to the above discussion. First, in 
the Wilson procedure of treatment the strong coupling domain of QCD one is to cut 
the perturbative integrals at the scale /i to ensure the reasonably small value of the 
coupling at /i and above. On the other hand, it is not necessary per se to take /i 
as low as possible even for purely perturbative calculations: the operator expansion 
can be applied to the perturbative fields as well, expressing the contribution of the 
low momentum region in the perturbative loop diagrams in terms of the perturbative 
matrix elements. Moreover, in such a way the contribution of this region is most 
simply evaluated and even computed (see, for example, papers @, where this 



idea was used in the heavy quark expansion); following this strategy one can easily 
go beyond the first loop, or with reasonable complication, even beyond the BLM 
approximation, which normally is difficult for complete analytical calculations. For 
example, the characteristic scale in 6 — > c is set by mi, — rUc or 2mc, which is rather 
large, and thus safely allows for treatment of the domain up to 1 GeV in this way. 
As an example, we can consider the expansion of the pole mass as the lowest 



eigenvalue of the rest frame Hamiltonian [0 obtained in l/mg approach to some 
order k in perturbation theory: 

<• - ^ I A„(.)-J^ ^ A(.) + ^ + ... (15) 

where the series runs in inverse powers of mQ(/i). According to Eq. (|15|), the per- 
turbative matrix elements A„(/i) are given in the BLM approximation by the corre- 
sponding coefficients of the expansion of J-'m{t) (it is defined in Eq. (|12[)) at small t: 
for 



^mif) — ^ Cut- 



one has 



n=0 

^ nil / t^^. { • , \ \ ^ 



TT Jo 



n+1 TT 
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For n = we have Ao(/i) = A(/i) and, therefore, its /i-dependence is indeed given 
by Eq. (^. Turning to the next term in l/mg, n = 1, we see that Tmif) does not 
contain the term ~ t^/^, i.e. Ci = 0. It means that the matrix element of the kinetic 
operator does not receive the perturbative contribution in the BLM approximation 
to any order in (foas/vr) when the cutoff is introduced in this particular way (this 
fact was pointed out in the somewhat different context in Ref. 0]). Literally the 



same technique can be applied to other quantities, not only well. 

The second comment concerns the concrete way of eliminating the low energy 
momentum region implicitly implied in our reasoning. Clearly, in general it does 
not reduces only to using the running quark masses: even after expressing the 
width in terms of mQ{fi) integration in Feynman graphs determining the width 
runs, strictly speaking, over all momenta starting fc^ = 0. Therefore formally an 
additional subtraction of the remaining contribution below /i^ is necessary. It is 
easy to see using the above arguments, however, that in this particular case the 
situation is rather simple: these extra terms can be neglected. Indeed, we can use 
the operator relation for the width [^, ^ 



Fsi oc ml{fi) ■ z 



pelt ({jTi-Jml)] fi 




2ml 



^pert 



"if. 



(17) 



where /x^ and /x^ are the expectation values of the kinetic and chromomagnetic op- 
erators normalized at point /i. The chromomagnetic operator does not mix with the 
unit one due to the different spin structure, and therefore its perturbative matrix 
element over heavy quark vanishes. As was mentioned just above, in the leading 
BLM approximation the perturbative value of //^ vanishes as well if one introduces 
the cutoff in the way adopted here. Thus one immediately concludes that the contri- 
bution of the infrared region to the perturbative coefficient of the leading operator 
in Eq. (p!?!), Cpert; is suppressed by at least the third power of jJ^/mq to any order 
in a<j in the BLM approximation 0. The BLM-type calculations of apert(/^) then 
very infrared stable. 

We see that the contribution of the domain below /z to the perturbative correc- 
tion factor expressed in terms of mQ(/i) is given by the corresponding perturbative 
contribution to the matrix elements of higher dimension operators only whose cor- 
rections to width scale like l/mq and higher powers of the inverse mass. On the 
other hand, nonperturbative corrections in the width are so far accounted for only 
through terms ~ l/rriQ. Thus as long as at the scale /i ~ 1 GeV the perturbative 
part of matrix elements is dominated by the condensate-like effects, the neglection 

■^This exactly corresponds to the result of Ref. which extended the statement of Ref. jij 
about the absence of terms ^ ^/mg through terms ~ I/uiq; however, the latter seems to be a 
feature of the BLM approximation only. 
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of the remaining perturbative corrections from the domain k"^ < fi^ is not only 
legitimate, but rather necessary. 

Let us note that in the exclusive B ^ D^*^ transition the dependence of the zero 
recoil amplitude on the masses enters only via the perturbative (and also nonpertur- 
bative) corrections. Strictly speaking, even here one should have used the running 
masses to calculate, say, the perturbative factors rjAy- Then, for example, to the 
second order in one would have 

VA^VA + c.m — —y{l-y)^c:^r]A-Cm(—) — ( 1 + ?/ + log ?/ ) (18) 
IT rric ay \ tt / rric \ ^ — y / 

y = rric/mb , fi <^mc . 

The effect, however, does not appear in the BLM approximation; numerically it 
is insignificant compared to nonperturbative corrections 0. The peculiarity of the 
determination of \Vqb\ from the total semileptonic width is that the latter in the tree 
approximation is proportional to the power of the quark masses, and therefore it 
does require the accurate treatment of the heavy quark masses. 



2 Implications for Determination of \Vcb\ and \Vub\ 

Let us briefly discuss the phenomenological implications. Eliminating the spuri- 
ous l/mQ renormalon which, in fact, is not present in the relation between differ- 
ent observables in QCD |^ by using the running quark masses instead of the pole 
masses, makes the perturbative corrections essentially smaller. Now we would get 
for rric/mh = 0.3 and m^, = 4.8 GeV 



^5 



m 



1927r^ \ TT 




rsi(& - c) ^ '^^^ [I - OM -^-'--^ + 1-7 ( - ) ) (19) 

rfifc = mf,(1.3 GeV) , rfic = "^£(1.3 GeV) . 

Addressing the CKM mixing parameters Vub and Vcb, one needs to consider the 
square root of the widths, and in particular the corresponding perturbative factors 



^From the theoretical perspective, however, this modification is profound: it represents the 
only linear in l/mg correction to rj^ as it is defined in HQET. To phrase it differently: unless one 



uses the running masses, i.e. introduces the corrections similar to those in Eq. (18), the known 
statements about the absence of corrections linear in 1 /mg to the zero recoil amplitudes from the 
low energy domain, are incorrect. (For the most general justification of the above conclusion see 
part 7.) The last term in Eq. ( |l8|) th us represents the counter-example to the Luke's theorem 
as it is formulated in HQET, see, e.g. |pO| ; this effect appears only beyond the BLM approximation, 
at the order and has not been noticed so far. 



9 




?7r = (r^i'^'^Vrsi) which are directly related to extraction of |V^fe|. In terms of m 
they then look as follows: 

r7r ^ I 1 - 0.31^^^^^ - 0.6 
V ^ 

r^r ^ (^1 - o.43 "-(V^^) . + o.8 ( ^ ) ) for 6 ^ c (20) 

These approximate expressions show the real sensitivity of extraction of \Vub\ and 
\Vcb\ to the perturbative effects. For example, the impact of the second order BLM 
correction is smallish. Moreover, because now the values of the heavy quark masses 
are essentially uncorrelated with a^, the first order coefficient indicates the actual 
uncertainty in \Vcb\ associated with the precise value of a^, which thus appears to be 
small and does not exceed the one occurring for the exclusive B ^ D* zero recoil 
amplitude. 

Using the most recent evaluation of we can now obtain the updated determi- 
nation of \Vch\ and the expression for \Vub\] although they are not really affected by 
the second order BLM-type corrections, we will use the calculated terms liter- 
ally to obtain the 'central' theoretical values. Two clarifying comments are in order 
beforehand. 

The above numerical values of the coefficients were shown for illustrative pur- 
poses only; in fact, the ratio mc/mb even for the one loop pole masses appears to be 
less than 0.3 when the l/mg expansion relations (see Eq. (^T]) below) are respected; 
this ratio decreases further for the running masses. In our numerical analysis we 
compute the values of the coefficients anew for proper values of rrib and rric/mb. 

Second comment concerns the value of rric. In principle, it can be accurately 
determined independently from charmonium sum rules, or from the semileptonic 
b ^ c spectrum itself At present, however, the most accurate estimate 

follows from the relation 

mb-mc = ^ — ( Z: Z~ ) + ^ I TT (^l) 




which has been used for accurate determination of |V"c6| in Refs. |T6|. This relation 
is formally valid for masses normalized at any legitimate point fi ^ Aqcd, and for 
our purposes it would be most convenient to take /z directly the same as used in 
Sect. 1, ~ 1 GeV. Strictly speaking, the result depends to some extent on the 
particular choice of /i due to terms ~ (a^/vr) /x^/mg (and higher in as/n and/or 
fi/mg) contributing to Eq. ( ^1]) in the form of the perturbative pieces of expectation 
values of higher dimension operators. In particular, one may be concerned with the 
terms ~ a<(/i^/m^ for large values of /i ~ rric. In fact, this effect is easily controlled 
to all orders in the framework of the BLM approximation. 



mt 
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mb{fi) - rricifi) + ^^^^ I dk\j^rn{^]-^m{^]] + 




(22) 



The variation of the mass difference does not exceed 20 MeV even when /z' varies 
up to 1.4 GeV. This dependence on the normahzation point is less than the error 
associated with the existing uncertainty in /x^ and can be discarded |. For the sake 
of definiteness we use the relation (PT|) at the scale 0.5 GeV; this literally decreases 



the value of — mc(/i) obtained without radiative corrections by only a few 

MeV. 

For our numerical estimates we use, as the central value, the strong coupling 
determined from the sum rules for bb production |^ as well: 

af^{l GeV) = 0^(2.3 GeV) = 0.336 ± 0.011 (23) 

and the value 

ml ~ mb(1.3 GeV) ~ m{;°^' - 0.560^(2.3 GeV) ■ 1 GeV = 4.639 GeV ; (24) 

in view of the numerical close proximity of the SV scale y^rnjfih to the above V 
scheme scale we expressed the perturbative corrections in terms of a^(2.3GeV) 
in the case of 6 ^ c decays irrespective of the exact values of the running quark 
masses, and for 6 — m decays use a^(4.8GeV); the second order coefficients are 
then adjusted appropriately. Again, this does not lead to any noticeable variation 
as compared to different choices. 

We then obtain for /x^ = 0.5 GeV^ using the normalization scale = 1 GeV : 



Kfel = 0.00458 



'Bt{B XJu)Y /1.6ps 



0.002 ) \ TB 
and 

IK.|= 0.0408 ( ' j (^) (25) 

It is important to emphasize that neither of the numbers above depend essentially 
on the exact value of the scale /z: varying it from (i.e. using the "two loop pole 
mass") to 1.35 GeV literally changes the value of \Vcb\ by only —1.2 percentage 
points. Using the running masses is only instructive in showing explicitly the small 
effect of higher order corrections; it is not mandatory in practice as long as one 
uses the same approximations in the theoretical expression for the widths and for 

®The suppression of corrections to the mass difference reflects the fact that they start with 
terms ~ (/i/mg)'^; this foUows from the absence of the perturbative contribution to /i^ in the BLM 
approximation. 
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determination of quark masses |T6[. The dependence is more pronounced in the 
case of 6 — >• M, but even here \Vub\ increases by only 2 percentage points when fj, 
is descended from 2 GeV down to 0.5 GeV. Let us finally note that discarding the 
second order perturbative terms ~ 6(0:5/77)^ altogether would literally lead to the 
explicit coefficients 0.00460 and 0.0410 in the above equations, respectively, provided 
one uses the same running masses normalized at 1 GeV. The change is less than a 
half of percentage point [|! 

Turning to the estimate of the actual theoretical accuracy of extracting \Vcb\ 
we note that the dominant uncertainty comes from the value of /i^: varying it by 
±0.1 GeV^ changes \Vcb\ by the factor (1 =1= 0.013) which, in fact, comes mainly from 
the related variation of m^, — m^. The dependence on the exact value of ml is rather 
weak: changing it by ±30MeV leads to the shift in \Vcb\ by =f0.6 percentage points. 

The main uncertainty in \Vcb\ is associated with the exact value of rub — rric. 
Changing it by ±30 MeV as compared to the one given by Eq. (^Tj) leads to the 
factor (1 =p 0.011). The relation (pl|) is exact to the order l/m^, however it may 
well be affected at the level of 20 MeV by terms ~ 1/m^. Therefore, at a percent 
level of accuracy one needs to estimate 1 / terms in the heavy quark expansion if 
relies on this mass relation for fixing iric. On the other hand, the value of — iric 
can be accurately determined from the semileptonic spectrum itself not appealing to 
the expansion in l/m^; then the nonperturbative corrections calculated 0] through 
terms l/rn^ are sufficient to determine the mass difference with the necessary ac- 



curacy p2[, and even when rric is not large. Combining the two methods one can, 
therefore, eliminate this source of uncertainty in \Vcb\ at a percent level in the both 
limits of large and small uic-, i.e. have the determination of \Vcb\ limited, on theo- 
retical side, only by terms oc 1/m^. 

Finally, the dependence on the exact value of is also rather moderate: changing 
a^^(lGeV) by ±0.02 emerges in the variation of \Vcb\ by ±0.6 percentage points 
only. 

Even better stability would hold for \Vub\'- the dependence on /i^ is very weak 
here, ±0.12 percentage points for every ±0.1 GeV^ in the latter. The similar change 
in ml by ±30 MeV generates the variation by the factor (1 =p 0.016) and the uncer- 
tainty of ±0.02 in ap(l GeV) translates into the factor (1 ± 0.009). 

The fit error bars quoted in Ref. are 

5ml = ±0.002 GeV , (5ap(l GeV) = ±0.011. (26) 

If used literally they would lead only to the negligible error in \Vcb\ less than 0.3%. 
It is clear, however, that for our purposes it cannot be taken at face value, even 
leaving aside the presently unknown value of yU^. There are other sources of (a^/vr)^ 

^Comparing with the previous estimate ||^, the central value of | V"c&| obtained there for the same 
input parameters Tsi{b —> c) and /ij as adopted in the present paper, reads \Vcb\ = 0.0405. In this 
respect it is worth clarifying that the numerical estimate in Ref. Eq.(25), was obtained using 
the V scheme coupling corresponding to the MS one quoted in jo); the latter was 15% larger than 
the central value adopted in the present paper. 
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perturbative corrections, both in the inclusive width and in the determination of 
ml in Eq. (^) which are probably more important. Of course, the above quoted 
error in is not relevant here as well, because ml cannot be viewed as the exact 
value of the running mass at a known scale. To be on the conservative side we 
feel necessary, before the dedicated analysis [^] for the running mass is completed, 
assign much larger uncertainty to mfe(lGeV), 6mi, ^ 50MeV (the uncertainty in 
this running mass in the analysis of sum rules for bb production seems to be about 
20 -T- 30MeV ^3[)' ^"^^ allow for the uncertainty in a2 for the inclusive widths of 
at least ±(3^5) as long as the complete two loop calculations for the width are 
not available. Adding these theoretical error bars one ends up with the current 
theoretical accuracy in \Vcb\ of about (3 3.5)% as a rather conservative estimate, 
provided /i^ is known. Anyway, it is quite possible that the theoretical precision 
cannot be reliably pushed below a percentage level due to potential preasymptotic 
corrections given by "exponential" terms limiting the applicability of duality for 
semileptonic decays of actual b hadrons ||24|| . 

The similar theoretical accuracy of the hypothetical determination of \Vub\ from 
the inclusive width r{B —>■ X„£z/) might seem to be extraordinary good, better than 
3%. In fact, rather sizable effects here may, in principle, come from terms ~ ^/m^, 
in particular due to generic Weak Annihilation processes possible in the KM 
suppressed semileptonic decays, which are expressed in terms of the expectation 
values of local four fermion operators PUI ^7\, B^, M ; their effect can be numerically 



enhanced here. To some extent it can be controlled by studying semileptonic KM 
suppressed decays separately for charged and neutral B mesons, and in particular 
in the end point region where the effect mainly originates from [^]. Therefore the 



determination of | Vub \ based only on the high energy part of the lepton spectrum can 
undergo strong higher order nonperturbative corrections that are not well known yet. 
If the total b —>■ u semileptonic width were accurately known then the uncertainty 
in \Vub\ would not exceed 5% level. 



3 Conclusions 

We have pointed out in this paper that the large (in particular for the 6 — > u channel) 
second order perturbative coefficients for the inclusive semileptonic widths of beauty 
particles are mainly associated with the similar contributions to the pole masses of 
heavy quarks when the widths are expressed in terms of The corrections be- 

come small if one uses theoretically well defined running masses normalized at the 
scale about 1 GeV; moreover, only these masses can and, as a matter of fact, are 
determined from experiment with necessary accuracy. In particular, the systematic 
shift in the values of \Vcb\ and \ Vub\ expressed in terms of the corresponding semilep- 
tonic widths, constitutes less than one percentage point when proceeding from the 
first order perturbative expressions to the ones where the corrections are calcu- 
lated using the BLM approximation, provided the value of mf,(l GeV) is fixed. Thus 
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the perturbative corrections, as well as nonperturbative effects, seem to be under 
good control at the level corresponding to a percent relative accuracy in \Vgb\. 

Even in the semileptonic decays of charm the consistent OPE treatment leads to 
the smallness of perturbative corrections: the large negative value of perturbative 
corrections found in Refs. in fact came from the momentum region near and 
below the infrared pole in the strong coupling. This contribution is therefore nu- 
merically irrelevant and is to be excluded from the calculations applied to charm; it 
is completely accounted for by nonperturbative effects within Wilson OPE. For this 
reason we cannot consider convincing the conjecture stated in Refs. 0] that the ap- 
parent numerical discrepancy of experimental semileptonic width of D mesons with 
theoretical expectations is associated with the uncontrollable nature of the pertur- 
bative series. The possibility to witness sizable violations of duality in this case 
seems to be more probable. 

The theoretical accuracy of calculations of total semileptonic widths of b particles 
appears to be very good when the input from the analysis of the bb threshold domain 
is used. One has, as the central values, at /i^ = 0.5 GeV^ 

= 0.00458 f5£(^^^ 
' ' I 0.002 




, , (Bi(B XJu) 

The main uncertainty in \Vcb\ at present comes from the exact value of /i^: 
/ ,(^^-0.5GeV2)\ / 6ml ^ 



iKbloc 1-0.013^^-2 — - — ^ . 1 - 0.006 . (28) 

' 0.1GeV2 ; \^ SOMeVy ^ ' 

The conservative estimate of other uncertainties associated, in particular, with not 
yet calculated part of the second order perturbative corrections, is about 3 percent- 
age points. The above expressions for \ Vcb\ practically do not differ from the previous 
estimate obtained in Ref. p. 

The theoretical accuracy of calculating Ts\{b — > u) is even better; there is no 
strong dependence on /i^ here, and the precision is better than 5% of the equivalent 
relative variation in \ Vub\- The uncertainty associated with the mass of b quark does 
not exceed 2 -h 3 percentage points and is not dominant. 

The analysis made in this paper suggests also that the impact of the third and 
higher order perturbative corrections, which can be computed straightforwardly 



within the BLM approximation |]lT], |T3|, |15|, is to be small and under good theo- 
retical control provided one treats the quark masses in the consistent way. Namely, 
if the pole mass for is used, its numerical value must be calculated in exactly the 
same approximation as applied to the calculation of the width. For example, one 
can impose the constraint that the one (or two) loop pole mass is not changed by the 
higher order corrections. More simple, and theoretically appropriate, approach is to 
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express the corrections in terms of the running mass at the scale around 1 GeV, which 
can be even better determined from experiment. (The usual technical problem with 
the precise definition of normalization point is absent in the BLM approximation.) 
The BLM corrections can then effectively sum the potentially dominant terms com- 
ing from high momenta by taking into account running of a^. Because the coupling 
in this region is relatively small, the running is not sharp and is well approximated 
already by the second term in the expansion of a^. For this reason one expects 
numerically small impact of the BLM corrections beyond blas/TT)"^ terms, and the 
obtained series must be under good numerical control. On the contrary, the effect of 
the low momentum physics is to be computed by means of the Wilson OPE rather 
than summing up perturbative series unstable in infrared. 

Note added: When this paper was in progress I was informed by V. Braun about 
computations of the effects of higher order perturbative corrections for semileptonic 



widths made within the technique of papers |jT3|, |T5|], Ref. p8|. The preliminary 
results reported agreed with the expectations that the actual impact of higher order 
corrections is smallish when the proper corrections in the quark masses are intro- 
duced. 
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b^c 


b^c 


b ^ u 


b ^ u 






0.2 


CLi 


0,2 


0.3 


-1.50 


-1.9 


-2.00 


-12.6 


0.5 


-1.38 


-0.9 


-1.72 


-9.5 


0.75 


-1.23 


0.1 


-1.38 


-6.3 


1.0 


-1.07 


0.9 


-1.03 


-3.7 


1.25 


-0.89 


1.6 


-0.70 


-1.5 


1.5 


-0.71 


2.1 


-0.35 


0.4 


2.0 


-0.32 


2.8 


0.32 


3.5 



Table 1: Dependence of the perturbative coefficients oi and 02 on the scale fi for 
rUb = 4.8 GeV and rric/mb = 0.3. The strong coupling is assumed to be defined 
in the scheme and normalized at ^JrajUb for 6 — * c and at for b ^ u . 



17 



